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c3 1. Introduction 



ABSTRACT 



Aims. We study stability of differentially rotating non-magnetic proto-neutron stars. 
Methods. Linear stability properties are considered taking into account neutrino transport. 

Results. If neutrino transport is efficient, the star can be subject to a diffusive instability that can occur even in the convectively stable 
region. The instability arises on the time-scale comparable to the time-scale of thermal diffusion. 

Conclusions. Hydrodynamic motions driven by the instability can lead to anisotropy in the neutrino flux since the instability is 
suppressed near the equator and rotation axis. 
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Newly born neutron stars are subject to various hydrody- 
namic and hydromagnetic instabilities shortly after core col- 
lapse. These instabilities can play an important role in enhanc- 
ing neutrino luminosities and increasing the energy deposition 
efficiency. Convection is probably among the best-studied in- 
stabilities in proto-neutron stars (PNSs) where it can be driven 
by both lepton and entropy gradients (see, e.g., Epstein 1979; 
Arnett 1987; Burrows & Lattimer 1988). The development of 
negative entropy and lepton gradients is common to many sim- 
ulations of core collapse supernovae (Bruenn & Mezzacappa 
1994; Bruenn, Mezzacappa & Dineva 1995; Thompson et al. 
2005) and evolutionary models of PNSs (Burrows & Lattimer 
1986; Keil & Janka 1995; Sumiyoshi et al. 1995; Pons et al. 
1999; Dessart et al. 2006). Numerical simulations also indicate 
the presence of instabilities in PNSs, however, conclusions re- 
garding the efficiency of turbulent transport and the duration 
of unstable phase depend on the code used and can differ sig- 
nificantly; compare, for instance, the simulations by Keil et al. 
(1996) and Mezzacappa et al. (1998). 

The nature of instabilities arising in PNSs has been consid- 
ered by a number of authors (Grossman et al. 1993; Bruenn & 
Dineva 1996; Miralles et al. 2000, 2002). The true stability cri- 
teria for non-magnetic and non-rotating PNSs (Miralles et al. 
2000) indicate the presence of two essentially different insta- 
bilities with the convectively unstable region surrounded by the 
neutron-finger unstable region, the latter typically involving a 
larger portion of the stellar material. Due to the cooling, the 
temperature and lepton gradients are progressively reduced, and 
both instabilities disappear completely after ~ 30 - 40 s. The 
growth time of convective and neutron-finger instabilities dif- 
fers substantially, and therefore the efficiency of turbulent trans- 
port in the convective and neutron-finger unstable zones can 
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be different as well. This difference is crucial for many MHD- 
phenomena in PNSs, such as turbulent dynamo action (see, e.g., 
Bonanno et al. 2003, 2005) and transport of the angular momen- 
tum. Note that fast rotation of PNSs can substantially modify 
convection and make convective motions anisotropic and con- 
strained to the polar regions (Fryer & Heger 2000, Miralles et al. 
2004). This mechanism is a natural way of creating anisotropic 
energy and momentum transport by convective motions, which 
only requires that the angular velocity be of the order of the 
Brunt-Vaisala or leptonic buoyant frequences. 

The hydrodynamic stability properties of proto-neutron 
stars, however, can be much more complex than is believed at 
the moment. Recently, for example, Bruenn et al. (2004) have 
argued that a new doubly diffusive instability can occur in an ex- 
tensive region below the neutrinosphere, which the authors refer 
to as "lepto-entropy fingers". This instability is driven mainly 
by the entropy equilibration due to a lepton fraction difference 
and lepton equilibration that is in turn caused by the entropy dif- 
ference and that may play an important role in enhancing the 
neutrino emission. 

The possible presence of the magnetic field and differential 
rotation in a core collapse supernovae favours a number of other 
instabilities that can also be important in enhancing anisotropic 
turbulent transport. For example, the magnetorotational insta- 
bility in the context of core collapse was first considered by 
Akiyama et al. (2003). The authors argue that the instability 
must inevitably occur in core collapse and that it has the ca- 
pacity to produce fields that are strong enough to affect, if not 
cause, the explosion. Thompson et al.(2005) constructed one- 
dimensional models, including rotation and magnetic fields, to 
study the mechanism of energy deposition. They explore several 
mechanisms for viscosity and argue that the turbulent viscosity 
caused by the magnetorotational instability is the most effective. 
The nonaxisymmetric magnetorotational instability in PNSs has 
been considered by Masada et al. (2006a). The authors have ob- 
tained the criterion and growth rate of instability and argue that 
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the nonhnear evolution may lead to enhancement of the neu- 
trino luminosity. The effect of neutrino radiation on this insta- 
bility was then considered by Masada et al. (2006b). Neutrino 
transport is rather fast in PNSs and can essentially modify the 
standard magnetorotational instability. 

In the present paper, we consider the effect of neutrino trans- 
port on the stabihty of differential rotation in PNSs, which can 
appear for a number of reasons. From theoretical modelling and 
simple analytic considerations, it is commonly accepted that the 
core collapse of a rotating progenitor leads to differential ro- 
tation for newly-born neutron stars (Zwerger & Mtiller 1997; 
Rampp et al. 1998; Liu 2002; Dimmelmeier et al. 2002; Miiller 
et al. 2004), mainly due to conservation of the angular momen- 
tum during collapse. It is also possible that PNSs are differen- 
tially rotating due to turbulent transport of the angular momen- 
tum during the convective phase. The boundaries of the convec- 
tive zone move inward when the PNS cools down and, in the 
region where convection stops, the angular velocity profile will 
evolve under the influence of shear viscosity and, what is more 
hkely, due to development of various instabilities caused by dif- 
ferential rotation. 

We show that this rotation can be unstable if it is accompa- 
nied by sufficiently rapid diffusion of heat and lepton number 
and if the angular velocity depends on the vertical coordinate. 
Contrary to the magnetorotational instability, the instability con- 
sidered in this paper does not require the presence of the mag- 
netic field. Turbulent motions caused by the instability can lead 
to anisotropic turbulent transport in convectively stable regions. 
The instability can be important when considering various dy- 
namo models, the transport of energy, and the angular momen- 
fiim in PNSs. 

The paper is organised as follows. In Section 2, we consider 
the basic equations governing the instability in differentially ro- 
tating PNSs and derive the dispersion equation. The stabiUty cri- 
terion and the growth rate of instability are treated in Section 3. 
A discussion of our results is presented in Section 4. 



2. Dispersion equation 

We now consider the PNS rotating with the angular velocity 
Q = Q.(s, z), where s, ip, and z are cylindrical coordinates. In 
the unperturbed state, the PNS is assumed to be in hydrostatic 
equilibrium. 



= G, G = g + Q}s 



(1) 



where g is the gravity and is the radial unit vector. We as- 
sume that the matter inside a PNS is in chemical equilibrium; 
thus, the pressure p is generally a function of the density p, tem- 
perature T, and lepton fraction Y - {tie + nv)/n, with tie and 
Hy the number densities of electrons and neutrinos, respectively, 
and n = np + n„is the number density of baryons. Therefore, we 
have Vln/7 = (dlnp/dlnp)TY(yinp + j3VlnT + 6VY), where 

and 6 are the coefficients of thermal and chemical expansion; 

= -idlnp/d\nT)pY,6 = -{din p/dY)pT . Then, taking the curl 
of Eq. (1), we obtain 



Gx{/3— +6VY\ = -e^s 



dz 



If O depends on z, hydrostatic equilibrium can be reached only if 
gradients of the temperature or lepton fraction have components 



perpendicular to G. If g » sQ^, then the variations of T and Y 
on a spherical surface r =const can be estimated as 



max 



1 



13 + 6 g dz 



Hydrostatic equilibrium can be satisfied if variations are small, 
or if sQ^ 5) 1 (we assume that 5Q/5z ~ Q/i). Typically, 
P+S ~ 0. 1 or larger in PNSs (see Miralles et al. 2002) and, hence, 
hydrostatic equilibrium can be reached for PNSs with dQ./dz + 
if the rotation period is longer than -1-2 ms. Most likely, the 
majority of PNSs satisfies this requirement and, therefore, the 
unperturbed state with £l{s, z) is justified for these PNSs. 

We consider axisymmetric short wavelength perturbations 
with the space-time dependence exp(yf - ik ■ r) where k = 
{ks,0,kj is the wavevector, \k ■ r\ » 1. Small perturbations 
will be indicated by subscript 1, whilst unperturbed quantities 
will have no subscipt. Since the growth time of the instability 
in PNSs is typically longer than the period of a sound wave 
with the same wavelength, we use the Boussinesq approxima- 
tion. Then, the linearized momentum and continuity equations 
governing the behaviour of small perturbations read as 



yv\ + 2£lxvi + e^sivi ■ VQ) 



*:vi = 0, 



p p 



(2) 
(3) 



where v\, p\, and p\ are perturbations of the velocity, pres- 
sure, and density, respectively, and v is the kinematic viscos- 
ity and the unit vector in the azimuthal direction. Since in 
the Boussinesq approximation, the perturbations of pressure are 
negligible, the perturbations of density can be expressed in terms 
of the perturbations of temperature and lepton fraction. 



Pi 



(4) 



During the initial evolution, the PNS is opaque to neutrinos, 
and the neutrino transport can be treated in the diffusion approx- 
imation (Imshenik & Nadezhin 1972). Then, the linearized ther- 
mal balance and diffusion equations read as 



Vi • — 

T T 



AVT ATi 

= kt^t +kyAYi, 



Fi + vi • vy : 



At^+AyAYi 



(5) 



(6) 



where At, Ay, kj, and ky are the corresponding kinetic coeffi- 
cients (see Miralles et al. 2000). We denote the superadiabatic 
temperature gradient as 



AVr = ( — I Vp-VT. 

\9PI.Y 



(7) 



The dispersion equation corresponding to Eqs. (2)-(6) is 
y'* + asy^ + a2y^ + a\y + ao = 0, (8) 
where 

a-i = 2a>y + a>T + Cl>y, 
a2 = ojtCOy - (OtyOJyt + 2ojyia)T + ojy) + cijl + q^ 

-(ojI + ojI), 



V. Urpin''^': Neutrino transport; and hydrodynamic stability of rotating proto-neutron stars 



3 



a\ = 2a»v(w2-wy - wjyWyr) + (o)t + ci)Y)(o)y + q ) 



-CJyicO^ + (Or) -ioAcOy 



-COyt 



In these expressions, we introduced the characteristic frequen- 
cies 



= yk', 



o>t = -fiA. 



(Oy = AyIc^, 
Oty = kyIc^, 

0)1 = SA- VY, 



where A = G - k(k ■ G)/k^. The quantities Wy, coj, and ojy are 
the inverse time scales of the dissipation of perturbations due 
to viscosity, thermal conductivity, and diffusivity, respectively; 
COyt characterises the rate of thermodiflftision, and coty describes 
the rate of heat transport due to the chemical inhomogeneity; 
ojg is the frequency of the buoyant wave; co^ characterises the 
dynamical time scale of the processes associated with the lepton 
gradient, and describes the influence of differential rotation 
on stability. 

If rotation is negligible (q^ « 0), then Eq. (8) reduces to a 

cubic one. 



7^ + i>2y^ + bij + ^0 = 0. 
where 

hi = COy + COT + COY, 
b\ = COtCOy — COtyCOyT + COyicOT + CiJy) - (Wg + Cijj), 



(9) 



bo = -CO^ I CjJy - -Uyt 



j - [cot - ^wryj 

+C0v(,C0tC0y ~ COty^Yt\ 



This equation is equivalent to the dispersion equation (17) de- 
rived by Miralles et al.(2000). 



3. Criteria and the growthi rate of instability 

Equation (8) describes four essentially different modes that exist 
in a differentially rotating PNS. The condition that one of the 
roots has a positive real part (unstable) is equivalent to requiring 
that at least one of the following inequalities be satisfied: 



^3 < 0, 
ao < 0, 

a^ai — a\ < 0, 
ai(fl3fl2 ~ ci\) ~ a^ao < 



(10) 

(11) 

(12) 

(13) 



(see, e.g., Aleksandrov et al. 1963). The quantity as is positive, 
therefore condition (10) will never apply. 
We have from condition (11) 



vco„ 



|/ly - |/lr j - vcol (kt - ^ky^ 
+{q^ + coI){ktAy - kyAt) < 0. 



(14) 



Note that, typically, the Prandtl number (the ratio of kinematic 
viscosity and thermal diffusivity) is small in PNSs (see, e.g., 
Thompson & Duncan 1993; Masada et al. 2006b), and the stabi- 
lizing influence of stratification is reduced. 

Substituting expressions for the coefficients ao,a\, ai, and aj, 
into conditions (12)-(13), we obtain rather cumbersome criteria 
of instability. However, these criteria can be much simphfied if 
we take into account that the frequencies have a different order 
of magnitude in PNSs. The dissipative frequencies COt, Cl>Yj Cl>tYj 
and COyt are approximately comparable to each other, and ol- 
ios' for perturbations with the wavelength of the order of the 
density lengthscale. The dynamical frequencies, cOg and co^, are 
typically much higher at ~ 10^ s '. Therefore, if analysing the 
conditions (12)-(13), we can restrict ourselves by terms of the 
highest order in dynamical frequencies and q^ since the angular 
velocity of PNSs is very likely much higher than 1 - 10 s~^ . The 
criteria of instability (12)-(13) read in the approximation as 



CO \v + Kt + 



co„\v + Av 



^/Ir j -col[v + Ay + ^/^yj 
+2vq^ < 0, 

-Col{v + KT- ^/Cyj 

+q^(KT + Ay) < 0. 



(15) 



(16) 



Conditions (15) and (16) represent neutron-finger and convective 
instability modified by rotation. If £2 = 0, these conditions agree 
with the criteria derived by Miralles et al. (2000). The effect of 
rotation is important for convective and neutron-finger instabil- 
ities if the angular velocity is comparable to the Brunt-Vaisala 
or leptonic buoyant frequences. In this paper, we concentrate on 
condition (14), which represents a new instability in PNSs asso- 
ciated to differential rotation and neutrino transport. 

We next consider the case of a moderately rapid rotation 
when the PNS is far from the rotational distortion and g » Q^i. 
Then, AVT, VF, and g are approximately radial, and 



co„ 



-p-liG.AVT) = -^co 



2 

'bv 



■ 6 -^iG■VY) = -^co■ 



io, 



(17) 



(18) 



where co^y = Pig ■ AVT) is the Brunt-Vaisala frequency, and 

'^io ~ ■ characteristic frequency of leptonic 

buoyant oscillations; ko is the 0-component of a wavevector, 
k^ — k^ cos^ 6 - 2ksky cos OsinO + k^ sin^ 0, is the polar an- 
gle. Substituting expressionns (17)-(18) into condition (14), we 
have 

k^ 

2 , 2 , e J. 2 
q + (^y + -j^ ^ <^b ^ 

where 

2 kt-/3ky/6 ^ v{Ay - SAtIP) 



(19) 



(Or, 



+ co; 



•w 



Ay ~ 6 At IP KtAy — KyAt 

Here, coi, is the characteristic buoyant frequency that is deter- 
mined by the Bruni-Vaisala frequency and the frequency of lep- 
tonic buoyant oscillations. 

Since ktAy - kyAt > from thermodynamics (see Pons et 
al. 1999), rotation-induced instabihty can occur in a stable strat- 
ification only if 



q 



d 



d 



K— -ks—\{s*Qr) <Q 



(20) 
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This condition is satisfactory if either the specific angular mo- 
mentum, s^il, decreases as the cyUndrical radius increases or Q 
depends on z such that dQ/dz ^ 0. 

Condition (19) depends on the direction of the wavevector 
and can be written as 



where 

A = 4QH iQj + ^(ol sin^ 6, 
B = sill+^ojlsinie, 
C = ^col cos^ e. 

In these expressions, we denote 



(21) 



(22) 



Since the dependence of F on the direction of k is quadratic, 
we can obtain that the relative maximum of F corresponds to 



1 ± 



(A-C)2 
' (A - C)2 + B2 



(23) 



The maximum value of F corresponding to these k^/k yields 
the following condition of instabiUty: 

2w^ + A + C ± Vfi^ + C)2 < 0. (24) 

The two conditions for instability then follow from the above 
expression: 



A + C + 2(ol = +^a>l+ IojI < 0, 



(25) 



or 



B'^ + {A-Cf>{A + C + lolf. 



where ^ = 4Q^ + is the epicycUc frequency. 

Inequality (25) can be fulfilled only if differential rotation 
satisfies the Rayleigh criterion, d{s^Q.)lds < 0. This criterion 
requires differential rotation that decreases very rapidly with the 
cylindric radius and usually is not fulfilled in stars. Condition 
(26) can be rewritten as 



- 4(A + lOy){C + iOy) > 0, 



or 



s^in^f + 2snl^iol sin 20 - Ak^^ojI - 
4colii^ + ^(4 + col) > 0. 



(27) 



(28) 



If stratification is stable (^w^ > 0) and rotation satisfies the 
Rayleigh stability criterion (k^ > 0), then the rotation-driven in- 
stability can occur only if dSl/dz + 0. Generally, such rotation 
can often be achieved in PNSs. Note, however, that even if the 
angular velocity depends on z, the instability does not arise near 
the rotation axis where s is small and near the equator where 
= 0. 

Since differential rotation is likely to be very strong in PNSs, 
we can estimate sO^ ~ Q^. The "viscous" frequency is negligi- 
ble in Eq.(28) if the wavelength of perturbations X = 2n/k satis- 
fies the condition 



(29) 



Since v ~ 10** - 10^ cm^/s in PNSs (see, e.g., Thompson & 
Duncan 1993; Masada et al. 2006), we can estimate Acr as lO'' 
cm, which allows a comfortable range of unstable wavelengths. 
For A > Acr, the necessary condition of instability is 

isiljf + 2s^1^ojI sin 26 - > 0. (30) 

If stratification is stable {^oj^ > 0), this inequaUty is equivalent 
to the conditions 



> ^ojI sin 20 + yjfojl sin^ 20 + Ak^^ojI, 
snl < ^ojI sin 20 - ■^^a)lsm^20 + 4K^^o)l. 



(31) 
(32) 



The first condition requires a large positive radial gradient of Sl^ 
that seems to be unlikely in PNS. On the contrary, the second 
condition can be satisfied if the vertical gradient of the angular 
velocity is negative. If rotation is sufficiently fast and /c^ ~ Q-^ » 



^w?, condition (32) yields 



Q » 2 y/^cob. 



(33) 



Since ^ ~ v/kj 1, condition (33) can be fulfilled in many 
PNSs even if Q < w/,. 

Many studies model rotation of the collapsing core by the 
angular velocity profile that depends on the spherical radius r 
alone, Q - fi(r) (see, e.g., Akiyama et al. 2003; Kotake et al. 
2004; Thompson et al. 2005; Sawai et al. 2005). Such a shellular 
rotation can be justified if the progenitor rotates with the an- 
gular velocity that depends on r alone (Monchmeyer & Mtiller 
1989). In this case, conservation of the angular momentum in 
the course of collapse leads straightforwardly to a shellular rota- 
tion of the collapsing core at the beginning of evolution, at least 
(see, e.g., Akiyama et al. 2003). For shellular rotation, we have 

= il^ sin 6 and = il^ cos where il^ = dD^/dr. Then, the 
criterion of instabiUty reads 



(26) 



sm20 ~ D.» 4^|^0Jb. 



(34) 



In this case, the instability occurs neither near the rotation axis 
nor near the equator but is strongly constrained to the polar angle 

~ 7T/4. 

The growth rate of instabiUty can be estimated from Eq. (8) 
if we assume that oji, ~ ciJbv > ^ and condition (30) is satisfied. 
Then, by making use of the perturbation procedure, we obtain 
for the unstable root, with accuracy in the lowest order in v and 



ao 

7~ 

ai 



(k: 



k'gColiAY - SAyIP) 



p 2 , 2 , 2 



(35) 



This root is positive (unstable) if criterion (19) is satisfied. Under 
conditions (29) and (33), expression (35) yields the following 
estimate for the growth rate of perturbations with kg ~ k: 



OJt — T. 



(36) 



Perturbations with small kg (that corresponds tok ^ kr) can grow 
even faster. In this case, we also should keep in ai the terms 
proportional to and obtain that y ~ coj. The growth rate of 
instability is rather high, particularly, for perturbations with a 
relatively short wavelength (but still satisfying condition (29)). 
Generally, y can be comparable to the growth rate of other diffu- 
sive instabilities in PNSs, such as the neutron-finger instability. 
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4. Discussion 

We have derived the criteria and growth rate of the rotation- 
driven instabihty in non-magnetic PNSs. This instabiHty is as- 
sociated with differential rotation and neutrino transport and can 
even arise in convectively stable regions. This instability is rep- 
resentative of a wide class of the so-called doubly diffusive in- 
stabilities known in stellar hydrodynamics (see, e.g., Acheson 
1978). In fact, this instability is the analogue of the well-known 
Goldreich-Schubert-Fricke instability for the case of neutrino 
transport (see Goldreich & Schubert 1968; Fricke 1969). The 
criterion of the Goldreich-Schubert-Fricke instability in ordinary 
stars is dQ./dz + 0. However, neutrino transport is qualitatively 
different from the radiative one even in diffusive approximation 
because, apart from transporting heat, it also involves transport 
of the lepton number (Eqs. (5) and (6)), which is important for 
the pressure balance. That is the reason the criterion of insta- 
bility in PNSs (30) is so different from the standard condition 
d^ldz + derived by Goldreich & Schubert (1968) and Fricke 
(1969). In PNSs, however, the instability also occurs only if the 
angular velocity depends on the vertical coordinate. 

During the early evolution stage, a rotation profile can be 
rather complex in PNSs and, therefore, the necessary condition 
of instability (33) can be fulfilled. For example, 2D calculations 
of the early rotation evolution of PNSs by Ott et al. (2006) indi- 
cate the presence of differential rotation in PNSs at 200 ms after 
bounce. Differential rotation is moderately strong in the inner 
region of a PNS with r < 10 km, where O varies by a factor 
1.5-3 depending on the model. In the outer region (r > 10 km), 
differential rotation is substantially stronger. 

The growth rate of the considered instability is small com- 
pared to that of convection; nevertheless, it is rather high and, 
in general, can be comparable to the growth rate of other dif- 
fusive instabilities such as the neutron-finger one. For example, 
even if rotation is slow and Q ~ 0.1 Wi, the growth time of per- 
turbations with a lengthscale close to the pressure scale-height 
is ~ Is. Perturbations with a shorter wavelength grow much 
faster. Most likely, then, the unstable motions reach saturation in 
PNSs. The turbulent velocity in saturation can be estimated by 
making use of the standard mixing-length approximation (see, 
e.g., Schwarzschild 1958). In the case of diffusive instabilities, 
this approximation agrees very well with numerical simulations 
(see Arlt & Urpin 2004). Then, the turbulent velocity with the 
lengthscale X is given by 

v(^) ~ Aj(X) ~ ^ ^, (37) 

where j(X) is the growth rate of perturbations with the wave- 
length A obtained from the linear theory of instability. If f2 ~ 
0.1-1 Wfo, then the velocity in the largest scale comparable to 
the pressure scale-height is ~ 10^ - 10^ cm/s. Such motions can 
essentially enhance, for example, the angular momentum trans- 
port that is important for the dynamo action in PNSs. The corre- 
sponding coefficient of turbulent diffusion is 

vr ~ — T, (38) 

and it can be comparable to the coefficient of heat diffusion that 
is largest among the diffisive coefficients in PNSs. Note that tur- 
bulent transport caused by the instability should be substantially 
anisotropic since it cannot arise in regions close to the equator 
and rotation axis. Such anisotropic transport can perhaps lead to 



an anisotropic enhancement in the neutrino flux during the initial 
second following the core-bounce. 

We have considered the instability of only short wavelength 
perturbations but, most likely, perturbations with a wavelength 
of the order of the pressure scale-heigh are unstable as well. 
Therefore, the instability can be the reason for a large-scale over- 
turn of the core region with a high energy density. This overturn 
can generally enhance the neutrino flux and make it anisotropic, 
which is crucial for the explosion mechanism since even a slight 
modification in the neutrino cooling/heating efficiency changes 
it qualitatively (Dessart et al. 2006). This is particularly impor- 
tant because the instability operates in the regions where convec- 
tion is suppressed and cannot influence the neutrino transport. 
For example, according to Dessart et al. (2006), convection at 
f < 1 s occurs in the region between 10 and 30 km, whereas the 
innermost region with the highest energy content is stable. The 
rotation-induced instability in the core could lead to motions that 
transport heat and lepton number to the outer region. The im- 
portant feature of the instability is that it produces large-scale 
departures from the spherical symmetry. This instability cannot 
arise near the equator and rotation axis and, therefore, turbulent 
transport is suppressed in those regions. 

Hydrodynamic motions induced by the instability can also 
contribute to the gravitational wave signals from the core- 
collapse supernovae. As recognised by Miiller et al. (2004), the 
dominant contribution to the gravitational wave signal is not 
produced by the core bounce itself but instead by the neutrino- 
driven convection in the postshock region and the Ledoux con- 
vection inside the deleptonizing proto-neutron star; note, how- 
ever, that Fryer et al. (2004) argue the opposite based on their 
results. Because of these contradictory results, the influence of 
convective motions on the gravitational wave signal seems to be 
an open issue. It is possible, however, that motions caused by 
the rotational instability, despite being slower than convective 
motions, can produce a significant gravitational signal, because 
the instability arises in the dense inner region of the PNS with 
the radius < 10 km where convection does not occur during the 
first second after the bounce. 
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